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The equations that govern the temporal evolution of two photons in the Schrodinger picture are derived, taking 
into account the effects of loss, group-velocity dispersion, temporal phase modulation, linear coupling among 
different optical modes, and four-wave mixing. Inspired by the formalism, we propose the concept of quantum 
temporal imaging, which uses dispersive elements and temporal phase modulators to manipulate the temporal 
correlation of two entangled photons. We also present the exact solution of a two-photon vector soliton, in order 
to demonstrate the ease of use and intuitiveness of the proposed formulation. 
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In quantum optics, the Heisenberg picture, where optical 
fields are treated as conjugate positions and momenta of quan- 
tized harmonic oscillators, is often prefeiTed, as it is easy to 
substitute the optical fields in classical electromagnetic prob- 
lems with non-commutative operators and obtain the Heisen- 
berg equations of motion. Once the operator equations are 
solved, one can then obtain various quantum properties of 
the optical fields via non-commutative algebra. However, the 
Heisenberg picture is not without shortcomings. It can be hard 
to analytically or numerically solve the complex or nonlinear 
operator equations without approximations. It is also diffi- 
cult to grasp any intuition about how the quantum correlations 
among the photons evolve until the Heisenberg equations are 
solved. These difficulties have led to a growing appreciation 
of the Schrodinger picture, where the photons are treated as 
an ensemble of bosons, and the evolution of the many-photon 
probability amplitude is studied. This arguably more intuitive 
approach has led to great success in the quantum theory of 
solitons lOl, where instead of solving the formidable nonlin- 
ear operator equations, one can obtain analytic solutions from 
the linear boson equations in the Schrodinger picture. The 
many-boson interpretation has been applied to the study of 
entangled photons as well, where the two-photon probabil- 
ity amplitude is shown to obey the Wolf equations by Saleh, 
Teich, and Sergienko (STS) |2]. Instead of treating the entan- 
glement properties of the photons and the optical propagation 
as two separate problems, with the STS equations, one can 
now use a single quantity, namely the two-photon amplitude, 
to keep track of the spatiotemporal entanglement evolution in 
free space. This is analogous to the Wolf equations, which re- 
formulate the laws of optics in terms of coherence propagation 
1. 

In this paper, we utilize the STS treatment of two pho- 
tons to study various temporal effects, in the hope that the 
Schrodinger picture would offer a more accessible interpreta- 
tion of temporal entanglement propagation for analytic or nu- 
merical studies of two-photon systems. Loss, group-velocity 
dispersion, temporal phase modulation, via an electro-optic 
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modulator for example, linear mode coupling, via a beam 
splitter or a fiber coupler for example, and four-wave mixing, 
in a coherently prepared atomic gas |4] for example, are all 
included in our proposed formalism, thus extending the STS 
model for use in many more topics in quantum optics, such as 
nonlocal dispersion cancellation |5, 6], fourth-order interfer- 
ometry |7], and two-photon nonlinear optics |4, 8]. The anal- 
ysis of a two-photon vector soliton, consisting of two photons 
in orthogonal polarizations under the cross-phase modulation 
effect, is presented in the final section, in order to demonstrate 
the ease of use and intuitiveness of the Schrodinger picture. 

Inspired by the formalism set forth, we propose the concept 
of quantum temporal imaging, which uses dispersive elements 
and temporal phase modulators to manipulate the temporal 
entanglement properties of two photons. Most significantly, 
we show that it is possible to convert positive time correla- 
tion to negative time correlation, or vice versa, using a tem- 
poral imaging system. This conversion technique should be 
immensely useful for applications that require negative time 
correlation, such as quantum-enhanced clock synchronization 
f9]. Although there have been theoretical 00,1301 and 
experimental 1 14] proposals of generating negative time corre- 
lation directly, they have various shortcomings compared with 
the conventional tried-and-true schemes that generate positive 
time coiTelation. Our proposed technique should therefore al- 
low more flexibility in choosing two-photon sources for quan- 
tum optics applications. 

The paper is structured as follows: Sec.|II|derives the equa- 
tions that describe the evolution of the two-photon amplitude 
in two separate modes. Sec. [Ill] introduces the principles of 
quantum temporal imaging, Sec. lIVI includes linear mode cou- 
pling in the formalism. Sec. IVl generalizes the formalism to 
two photons in more than two modes. Sec. IVII includes the 
effect of four- wave mixing, and Sec. I VIII presents the exact 
solution of a two-photon vector soliton. 



II. TWO PHOTONS IN TWO SEPARATE MODES 

Let us first consider two photons in two optical modes, 
such as two polarizations, two propagation directions or two 
waveguide modes. The corresponding two-photon wavefunc- 
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tion IS 



|'P)=Ci2|l,l)+Ci,|2,0)+C22|0,2) 



(1) 



where the constants Cj^'s are the overall amplitudes of the 
quantum states, |1, 1) is the quantum state in which one pho- 
ton is in each mode, |2,0) is the state in which both photons 
are in mode 1, and |0, 2) is the state which both photons are in 
mode 2. The positive-frequency forward-propagating compo- 
nent of the electric field in each mode is given by HJll^ 



(■ / dco (- 

Jo 



fi(or]j{co) 



yAn£Qc[nj{co)YSj 
dj{z,C0)exp{-im), (2) 

where «/ is the complex, frequency-dependent refractive in- 
dex in mode j, rjj is the real part of nj, S is an area of quan- 
tization in the x-y plane, and dj is the photon annihilation op- 
erator, related to the corresponding creation operator via the 
equal-space commutator fTstildl . 

[flXz,«),4(z,a)')]-5(«-«'), ; = 1,2. (3) 

In the Heisenberg picture, the creation and annihilation oper- 
ators evolve according to the following equations L15..16.1 . 
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(5) 



where pij is the imaginary part of nj, and fj is the Langevin 
noise operator, satisfying the commutation relation. 



[fj{z.co)J]{z!,co')]^8{z~z!)5{co~co'). 



(6) 



To proceed, we replace conj(co)/c by the following phe- 
nomenological approximation L17i] . 
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where aj = 2Im[kj(coo)] is the loss coefficient, j3„; = 
d"Re[kj{co)]/dco"\(o=(oo is the «th-order dispersion coeffi- 
cient, and Arij encompasses any other refractive index per- 
turbation. Defining the slowly-varying envelope operators as 

Aj{z,t) ^exp{-ij3ojZ + icoot) [ dj{z,co)exp{~icot), 

Jq v27t 

(8) 



where COq is the carrier frequency of the two modes, one can 
obtain two evolution equations for the envelope operators. 



y^Ai{z,t) = iKi(t,i^'^Ai{z,t) 



dz' 
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where Fj is defined as, 

Fjiz.t) ^exp{~ipojZ + ia)ot)x 

dco ,/2C0fij{C0)\7 



V27t V c 



'fj{z,co)exp{-icot), (12) 



and Kj is the complex wavenumber for the slowly-varying en- 
velope. Arij can explicity depend on time, if the perturbation 
is much slower than the optical-frequency oscillation so that 
an adiabatic approximation can be made, such as in an electro- 
optic modulator. 

We now define the two-photon probabiUty amplitudes as 



Wi2{z,t,z!,t')^{Q\A:{z,t)A2{z!,t')\^), 



(13) 



V/ii(z,f,z',f') = -^(0|Ai(z,f)Ai(z',OW, (14) 
W22{z,t,z',t') = -^(0|A2(z,f)A2(z',f')l'I'>- (15) 

The physical significance of each amplitude y/j^ is that its 
magnitude squared gives the probability density, Pji^, of co- 
incidentally measuring one photon in mode /' at (z,f) and an- 
other photon in mode k at {z! ,t'). 



P,-i(z,f,z',f') = |r,^(z,f,z',f')l^- 



(16) 



Temporal entanglement is defined as the irreducibility of 
|V/i2p into a product of one-photon amplitudes in the form 
of a{t)b{t'). This means that the probability of detecting a 
photon in mode 1 at time t is correlated to the probability of 
detecting a photon in mode 2 at f'. The most popular ways 
of generating entangled photons are spontaneous parametric 
down conversion 1 18] and four-wave mixing 1 191 . where the 
wave mixing geometry and the spatiotemporal profile of the 
pump beam determine the initial ^\2- 

To obtain the evolution equations for the two-photon am- 
plitude \j/i2{z,t,z' ,t') in the Schrodinger picture, we employ 
the same trick as in Ref. yj. First we multiply Eq. with 
A2(z',f') and Eq. J10> with Ai(z,f) to produce two equations, 

^Aiizj)A2{z\t')^iKi(^tJ^yiiz,t)A2{z',t')+FiA2, 

(17) 

^A 1 (z, f )A2 {z',t') - iK2 (f', A 1 (z, f )A2 iz' ,t')+F2M. 

(18) 

Using the definition of 1/12 in Eq. ( I13> and assuming that the 
thermal reservoirs are in the vacuum state so that the Langevin 
operators evaluate to zero when applied to the wavefunction 
I20ll . a pair of equations in terms of 1/12 are derived, 

^VAi2(z,f,z',f')-/^:i(f,/^)v^i2(z,f,z',f'), (19) 

^ri2(z,f,z',f') = '^^2(f',/^)ri2(z,f,z',f')- (20) 

Equations ( 119) and (|20} are the temporal version of the STS 
equations |j2J, including the effects of loss, dispersion and 



3 



phase modulation. They can also be written in the frequency 
domain as 



(j)niz,Q.,z' = J dtj dt' \i/niz,t,z',t')x 

exp{iQ.t + iQ.'t'), (21) 
0niz,^,z',a')^iKi(^-^,ayn{z,^,z',a'), (22) 
1 d 



dz 



— 0i2(z,n,z',n') = /^:2(-7^,ii'j(/'i2(z,n,z',n'). (23) 

For entangled photons, because 1/12 or 0i2 cannot be separated 
into a product of one-photon amplitudes, distortions experi- 
enced in one arm can coherently add to the distortions expe- 
rienced in the other arm, leading to various nonlocal quantum 
effects. 

For example, consider group-velocity dispersion only, the 
output 012 is given by 



0i2(z,n,z',n') 

= exp (i^n^z + iPn^'z' + '-^Slh + '-^a'h') x 
(?)i2(0,n,0,n'). (24) 



ip2 



If the photons are initially entangled with negative fre- 
quency correlation, 0i2(O,i2,O,i2') can be approximated by 
(p{Q.)5{Q. + Of). Ignoring the unimportant linear spectral 
phase, the output is 

(;)i2(z,n,z',n')=exp[— (j32iz + j322z')j0(ii)5(i2 + a'). 

(25) 

Hence if j32iz — —jhiz', the dispersion effects in both arms 
can nonlocally cancel each other, as originally discovered by 
Franson 



III. QUANTUM TEMPORAL IMAGING 

In the Schrodinger picture, the two-photon amplitude 
evolves under temporal effects. Since the entanglement prop- 
erties of the photons are contained in the two-photon ampli- 
tude, the Schrodinger picture allows one to use the temporal 
effects to engineer the entanglement. 

First, consider the evolution of the two-photon amplitude 
when one of the modes, say mode 1, is subject to group- 
velocity dispersion. 



(26) 



dz dt 2 dfi ' 

xifnmt,z',t') = dTbiit-T)xi/n{0,T,z\t'), (27) 

-/(f-j3nL-T)2^ 



exp 



2j32iL 



(28) 



Group-velocity dispersion is well known to be analogous to 
Fresnel diffraction. 



Next, consider a quadratic temporal modulation of refrac- 
tive index imposed on mode 1 by a short or traveling-wave 
electro-optic modulator. 



d\jfi2 ikoAn2{t — to)^ 



dz 

Vl2{l,t,z',t') 



Vn, 

q{t)\t/n{0,t,z\t'), 



q{t) = exp 



ikQAn2l I ^ ^2 



(29) 
(30) 
(31) 



Quadratic temporal phase modulation is analogous to a lens. 
A«2 is assumed to be a constant, and Iq is the time delay of the 
modulation. Kerr effect by a co-propagating classical pulse 
would also suffice. 

Two dispersive elements and a quadratic phase modulator 
in-between form a temporal imaging system, which has been 
well studied in the classical domain f2Tll . Suppose that the 
photon in mode 1 propagates through the first dispersive ele- 
ment, with an effective dispersion coefficient ^21 and effective 
length L, then passes through a time lens with refractive index 
modulation An2{t — ?o)^/2, and finally propagates through the 
second dispersive element, with an effective dispersion coeffi- 
cient 1^21 and effective length L'. The output two-photon am- 
plitude can be expressed in terms of the input as 

Vfl2{z,t,z',t')= f dx f dx' b\{t-T')q{z')bi{z'~z)x 

(32) 



V/l2(0,T,z',f') 

.27r^2i 



-i{t-j3[,L'-T'f 



When the "lens law" for the time domain is satisfied, 

= knArioL 



P21L j3^iL' 
the impulse response of the system becomes 

/j(f,T)= / dV b[{t-T')q{T:')bi{T' -t) 



(33) 



(34) 



(35) 



exp 



/(T + j3nL)2 



^i{t~ii[,L'f 



2J3'L' 



exp 
exp 



2j32iL 
A 13^, L' ^ j32iL 



dz'P{^)x 

in 



koAn2lto 1 t' 



■'21^ 



(36) 



where P(T'/r„) is the normalized temporal aperture function 
of the time lens that can be used to describe any deviation of 
the actual temporal phase modulation from the ideal quadratic 
profile, such as truncation or higher-order phase modulation, 
and Ta is the aperture width. If 



Ta » 



(37) 



where Tq is the smallest feature size of 1/12 along the t axis, the 
integral in Eq. i36\ can be approximated by a delta function. 
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We then arrive at the input-output relation for the two-photon 
amphtude, 



M 



M 



,z',t'), 



td = P'nL'+MPnL+il-M)to, 



(38) 
(39) 
(40) 



where an unimportant quadratic phase factor is omitted, 
is the time delay of the system, and M is the magnification, 
which can be positive or negative depending on the signs of 
Pi and j3^. 



Positive-Time-Correlated T Electro-Optic J^' 
Two-Photon Source ^ ^ Modulator , 



Pump 
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FIG. 2: (color online). A quantum temporal imaging system for 
quantum-enhanced clock synchronization. 




FIG. 1: (color online). Two-dimensional sketches of the two-photon 
probability amplitude before and after one of the photons is time- 
reversed. Uncertainty in arrival time difference is transformed to un- 
certainty in mean arrival time. 

The most interesting case is when M = — 1, and one of the 
photons is time-reversed. If the two photons are initially en- 
tangled with positive time correlation, v/12 can be written as 

xi/n{0,t,z',t')^a{t)b{t-t'), (41) 

where b is assumed to be much sharper than a. After photon 
1 has passed through the temporal imaging system with M — 
-1, 

Wn{z,t,7:,t')^a{ta-t)b(t^-t-t'). (42) 

The photons hence become anti-correlated in time. See Fig.Q 
for an illustration of this process. Since most conventional 
two-photon sources generate positive time correlation, but 
negative time correlation is desirable for many applications, 
one can use the temporal imaging system to convert the for- 
mer to the latter In particular, using the aforementioned tech- 
nique for the specific application of clock synchronization, the 
sub-classical uncertainty of arrival time difference, (f — f')/2, 
can be converted to a sub-classical uncertainty of mean arrival 
time, [t +f')/2, leading to a quantum enhancement of clock 
synchronization accuracy by a factor of \/2 over the classi- 
cal limit. In practice, the clock can be synchronized with the 
electro-optic modulator, so that the mean aiTival time is con- 
trolled by fo and thus the clock. The proposed setup is drawn 
in Fig.|2l 

The fidelity of time reversal is limited by parasitic effects, 
such as higher-order dispersion and phase modulation, and 



the temporal aperture r„, which adds a factor ~ jSji^'/^a 
the width of 1/ along the t axis and increases the overall un- 
certainty of the mean arrival time. The ultimate limit, apart 
from instrumental ones, is set by the failure of the slowly- 
varying envelope approximation, which only concerns ultra- 
short pulses with few optical cycles. 

Besides the above application, one can also convert neg- 
ative time coiTelation, which can be generated by ultrashort 
pulses for improved efficiency fill fl3l l23l . to positive time 
correlation. As evident from Eq. J38> . any desired correla- 
tion can actually be imposed on already entangled photons, 
by multiplying the original correlation with a factor of 1 /M. 

As group-velocity dispersion and temporal phase modula- 
tion play analogous roles in the time domain to diffraction 
and lenses, one can use Fourier optics |23], temporal imaging 
HHl . and quantum imaging |24] techniques to design more 
complex quantum temporal imaging systems. 



IV. TWO PHOTONS IN TWO LINEARLY-COUPLED 
MODES 

Suppose that the two modes are now coupled to each other, 
via, for example, a beam splitter or a fiber coupler Equations 
(|9j and ( I10> become coupled-mode equations, 

(^-//i:i)Ai =/Jc(z)A2+A, (43) 
- iK'^)A'^ - iK* (z')A'i + Fi (44) 

where K is the coupling coefficient, and for simplicity the cou- 
pling is assumed to be co-directional. The primes denote the 
evaluations of the functions at {z' ,t'). Any phase mismatch 
can be incorporated into K as a z-dependent phase. 

Procedures similar to those in SeclUproduce four coupled 
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equations for xf/n, \j/22, and xj/u, 

(^^^ - iKi'^V2\iru{z,t,z' ,t') ^ iKxifniz' ,t' ,z,t), (45) 

(^-^-iK^^V2w22{z,t,z'y)^iK*'wi2{z/,z,t), (46) 

- iKi^ wn{z,t,z!,t') = iKV2\f/22iz,t,z!,t'), (47) 

(^-^-iK!^^Wn{z,t,z',t')^iK*'V2xifniz,t,z\t'). (48) 

Any pair of Eqs. i46\ and (147 1 or Eqs. ( I45> and ( I48t can be 

combined to yield a single equation for y/^. 



-JC(z)jC*(z')ri2(z',f',Z,f) 



(49) 



Equation iA9i allows one to calculate the coupled-mode prop- 
agation of two photons in terms of 1/12 only, given the initial 
conditions of v/12, Yn, and v/22- Wii ™d \if22 can then be 
obtained from Eqs. J47t and ( I48> after 1/12 is calculated. 

To obtain some insight into Eq. ( I49> . consider only constant 
mode coupling, so that Eq. ( I49> becomes 



^^V^12(z,f,z',f') = -K-Vl2(z',f',Z,0- 



(50) 



The solution is 

V/I2(z,f,z',f' 



cos(k-z) cos( fcz') 1/12(0, f,0,f')- 
sin(K-z) sin(K-z')v/i2(0,f',0,f) + 

ism{Kz) cos(K-z')\/2v/22(0,f,0,f')+ 

/cos(K-z)sin(K-z')V2v/ii(0,f,0,f')- (51) 

At the coupler output, z~z! —L, 

\l/niL,t,Ly)^T\irniO,t,0,t')-R\ifn{0,t',0,t)+ 
iVTrR\i/22{0,t,0/)+ 
iV2JR\i/n{0,t,0,t'). (52) 

where T ~ cos^(fcL) and R = I — T = sin^(fcL). If we have 
one photon in each mode initially, only the initial condition of 
y/i2 is non-zero, and 

\\fn{L,t,L,t') ^ T\i/i2{0,t,Q,t') -/?v/i2(0,f',0,f). (53) 

From Eq. i53\ . one can see that the output amplitude is the 
destructive interference between the original amplitude and 
its replica but with the two photons exchanging their positions 
in time. In particular, for a 50%-50% coupler, T ^ R — I /2, 
complete destructive interference is produced if the two in- 
put photons are temporally indistinguishable. See Fig.^for a 
graphical illustration of the destruction interference. The in- 
troduction of variable distinguishability to photons, in order 
to produce varying degrees of destructive interference of 1/12 
via a beam splitter and to measure the two-photon coherence 
time, is the basic principle of the Hong-Ou-Mandel interfer- 
ometer I^TJ. 



'/'„(o,7',o,o 




FIG. 3: (color online). The quantum destructive interference via a 
coupler is determined by the overlap (dark grey area) of the two- 
photon amplitude v/i2(0,f,0,f') with its mirror image with respect to 
the / + /' axis, v/i2(0,f',0,f)- 



V. TWO PHOTONS IN MANY MODES 

If the two photons are optically coupled to more than two 
modes, such as four modes for two polarizations in each of the 
two propagation directions, or modes in an array of fibers 
coupled to each other, one in general needs A^(A^+ l)/2 two- 
photon amplitudes to describe the system. The propagation of 
the amplitudes in many modes is described by the following, 

d 



\ + 8jkWjk{z,t,z!,t' 



i L ^ji V^ + SikVik {z, t,z',t' 



where 



¥jkiz,t,z\t') = \l/kj{z',t',Z,t), Kjk 



kj- 



(54) 



(55) 



Further simplications can also be made if any of the coupling 
terms is zero. 

For example, let there be four modes; mode 1 corresponds 
to arm 1 with x polarization, mode 2 corresponds to arm 2 with 
X polarization, mode 3 corresponds to arm 1 with y polariza- 
tion, and mode 4 corresponds to arm 2 with y polarization. If 
only the same polarizations are coupled, the two-photon equa- 
tions are 



-iKx ^ 








\ 





\ 

( Vly/u Wn 

1/21 \/2V/22 



d 

Tz 

-iKy 



-iKy 



dz 
ri3 
W23 



W2-X 

1/31 ^32 \/2l/33 ^34 

V l/tl '/42 </43 V21/44 ) 



(56) 
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The following solution for the orthogonally polarized ampli- 
tudes can be obtained. 



/ V^i3(L,f,L,f') 
W24{L,t,L,t') 
\lfu{L,t ,L,t') 

\ VA23(L,f,L,f') 



\/TxTy — ^ RxR y 
— ^ RxR\ \/TxTy 

V/i3(0,r,0,O 

V/24(0,f,0,f') 

V/i4(0,f,0,f') 
r23(0,f,0,f') 




y/R^ 



(57) 



where 7^.,. ~ cos^(jQjL) and Rx.y = 1 — r,.y. In particular, if 
only the initial condition of 1/14 is non-zero. 



Wn{L,t,L,t') = /y7^v'i4(0,f,0,f'), 
y/24{L,t,L/) = /y7^i/i4(0,f,0,f'), 
y/i4{L,t,L,t') = JT^.yfu{0,tAt'), 



\l/23{L,t,L,t') = ^y^R^.\l/i4{0,t,0,t'). 



(58) 
(59) 
(60) 
(61) 



The singlet state for orthogonally polarized photons is pro- 
duced ifr, = r,, = 1/2 m. 



VI. FOUR- WAVE MIXING 



As envisioned by Lukin et ai, the third-order nonlinear ef- 
fects among two photons can become significant in a coher- 
ently prepared atomic gas |4]. The coupled-mode equations 
(143 > and (I44> then become nonlinear. 



(62) 



(63) 



— -iKijAi= iKA2 + /yAjAiAi + /T7A2A2A1 + 



iXAiAiAl+Fu 
- iK^)A'2 = iK*'A[ + iYA'lA'2A'2 + /77A'Ja'iA^- 



iX*'A[A\A'l+Fi, 



where 7 is the self-phase modulation coefficient, rj is the 
cross-phase modulation coefficient, and x is the four-wave 
mixing coefficient. If we define equal-space two-photon am- 
plitudes as the following. 



three linear coupled-mode equations for the two-photon am- 
plitudes can be derived, 

( ^ - iKi - iK[ ) V2 v/i i=iK\j/2i+ iyS (f - 1') V2\j/i 1 + 

ix5it^t')V2w22, (65) 
- iK2 - iKl2^ V2V/22 = iK*W2\ + 178(1 - t')V2\ir22+ 
/X*5(r-f')V2v/n, (66) 
- iKi - iK^'j 1/12 = /K-*\/2v/ii + /K-\/2v/22+ 

/T75(f-f')ri2. (67) 



The advantage of the Schrodinger picture is most evident here; 
whereas in the Heisenberg picture one needs to solve nonlin- 
ear coupled-mode operator equations such as Eqs. ( I62t and 
j63> . in the Schrodinger picture, one only needs to solve lin- 
ear equations such as Eqs. i65\ to ( I67t . which are similar to 
the configuration-space model applied to the quantum theory 
of soHtons 1TII26I1. 

The delta function 5{t — f') couples the two subspaces of 
Wi2{z,t,t'), SO entanglement can emerge from unentangled 
photons f^. To see this effect, assume that we only have four- 
wave mixing, so that Eq. (I67> becomes 

^\l/i2iz,t,t')^iri5it-t')\(fi2iz,t,t'), (68) 
oz 

which yields 

V/i2(L,f,f') =exp[/T7L5(f-f')]ri2(0,f,f')- (69) 

If the nonlinearity has a finite bandwidth Aco, the delta func- 
tion in time should be replaced by a finite-bandwidth function, 
for example a sine function. 



VAi2(i,f,f') =exp| 



irjL 



n{t-t') 

V/I2(0,f,f'). 



Act) 



(70) 



(64) 



Eq. ( I70> is the exact solution of the two-photon amplitude 
under the cross-phase modulation effect, while Eq. (7) in 
Ref. |4], presumably derived in the Heisenberg picture, is only 
correct in the first-order As \if\2{I^,t,t') cwnolht written as a 
product of one-photon amplitudes even if the two photons are 
initially unentangled, entanglement is generated. The phys- 
ical interpretation is that the two input photons act as pump 
photons to the spontaneous four-wave mixing process and are 
annihilated to generate two new entangled photons. 

Unlike temporal imaging techniques, which can only ma- 
nipulate the two-photon amplitude along the horizontal axis 
t or the vertical axis cross-phase modulation allows some 
manipulation of the two-photon amplitude along the diagonal 
time-difference axis, t — t' . Unfortunately, cross-phase modu- 
lation by itself cannot generate any temporal correlation, as it 
only imposes a phase on the two-photon temporal amplitude. 
In order to have more control along the f — f ' axis, one can 
combine the effects of cross-phase modulation and dispersion, 
as shown in the following section. 
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VII. TWO-PHOTON VECTOR SOLITONS 

In this section we study a toy example, namely, a soliton 
formed by two photons in orthogonal polarizations exerting 
cross-phase modulation on each other 1 17|. Although similar 
studies of two photons in the same mode under the self -phase 
modulation effect have been performed in Refs. |8], cross- 
phase modulation offers the distinct possibility of entangling 
two photons in different modes. 

Consider the case in which two polarizations have the same 
group-velocity dispersion, so that ^21 — P22 — p2, and there 
is one photon in each polarization. The evolution equation for 
\lfniz,t,t') is 



f- 



d_ 
dt 



ip2(d^ 



Defining time coordinates in a moving frame, 

we obtain the following equation for v/i2(z, T, t'), 

d d d \ 
3-+A— -A— 
dz ax dx' I 

ip2{d^,^ 



Wn- (71) 



(72) 
(73) 



-/tj5(t- t') 



(74) 



Equation MA\ is a simple linear Schrodinger equation, de- 
scribing a two-dimensional "wavefunction" 1/12(2, T,t') in a 
moving frame subject to a delta potential. To solve for 
explicity, we define new time coordinates. 



x+ = 



T_ = 



X-X' 



(75) 



Eq. J74t then becomes 



•5 . d 



ip2 
4 



\dxl 



dxl 



IT] 



5{x- 



(76) 



As evident from Eq. M6\ . the cross-phase modulation effect 
only offers confinement of 1/12 along the time difference (T_) 
axis, but not the mean arrival time (T+) axis. 
The only bound-state solution of y/^ is 



ri2(z,T+,T_) 

= exp 



A^ 

J2 
2A 



exp - I + i-^^- ) 

da 



^ (£2) exp - iQ.x+ + ^ii^z) • (77) 



The delta potential enforces S to take on the following value, 

(78) 



where 77 and P2 must have opposite signs. The final solution 
of v/12 in the frame of x and x' is therefore 

Wn{z,x,x') 



- exp 
exp 



./tj2/4 + a2\ ■ 





1 




2p2 



t-t'| + /-^(t-t') 

P2 



.2n 



(n) exp 



(79) 



The two-photon coherence time of a vector soliton is fixed, but 
the average arrival time is still subject to dispersive spreading 
and becomes increasingly uncertain as they propagate. See 
Fig Elfoi" an illustration. Hence, a two-photon vector soliton 
generates temporal entanglement with positive time correla- 
tion as it propagates. Similar to the idea of soliton momen- 
tum squeezing |27], one can also adiabatically change 77 or P2 
along the propagation axis to control independently the two- 
photon coherence time. 

Notice that the center frequencies of the two photons are 
shifted slightly, by an amount of ±A/j32, to compensate for 
their group-velocity mismatch, so that they can co-propagate 
at the average group velocity. This is commonly known as 
soliton trapping Ill7ll . 



FIG. 4: (color online). Quantum dispersive spreading of mean ar- 
rival time of a two-photon vector soliton. The cross-phase modula- 
tion effect only preserves the two-photon coherence time, giving rise 
to temporal entanglement with positive time correlation. One can 
also manipulate the coherence time independently by adiabatically 
changing the nonlinear coefficient along the propagation axis. 

If the nonlinearity has a finite bandwidth, then the potential 
becomes a finite-bandwidth function like the one in Eq. MQl . 
and multiple bound-state solutions can be obtained via con- 
ventional techniques of solving the linear Schrodinger equa- 
tion. 
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Vni. CONCLUSION 

We have derived the general equations that govern the tem- 
poral evolution of two-photon probability amplitudes in dif- 
ferent coupled optical modes. The formalism inspires the con- 
cept of quantum temporal imaging, which can manipulate the 
temporal entanglement of photons via conventional imaging 
techniques. The theory also offers an intuitive interpretation 
of two-photon entanglement evolution, as demonstrated by the 
exact solution of a two-photon vector soUton. To conclude, we 



expect the proposed formalism to be useful for many quantum 
signal processing and communication applications. 
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